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Abstract 

Eguchi-Hanson solitons are odd-dimensional generalizations of the four-dimensional 
Eguchi-Hanson metric and are asymptotic to AdSs/Z p when the cosmological con- 
stant is either positive or negative. We find soliton solutions to Lovelock gravity in 5 
dimensions that are generalizations of these objects. 
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1 Introduction 



The AdS/CFT correspondence has played a pivotal role in modern theoretical physics for 
more than a decade because of its conceptual importance and its broad applications pQ. It 
conjectures that for some gravitational theory on a D-dimensional spacetime of the form 
AdSd x Vl-D-d where M is a compact space of D — d dimensions (D = 10 for string-theory 
and D = 11 for M-theory) then there also exists a mathematically equivalent (d — 1)- 
dimensional Conformal Field Theory (CFT). As a consequence, solutions to the vacuum 
Einstein equations in (d — l)-dimensions have particular importance since they furnish both 
ground states for string theory and background manifolds for the dual CFT. 

An interesting case in point are generalizations of the Eguchi-Hanson (EH) metric. EH 
metrics are solutions of the four-dimensional vacuum Euclidean Einstein equations, and 
can be regarded as a special case of the Atiyah-Hitchen metric. This latter metric can be 
embedded in M-theory to furnish new M2- and M5-brane solutions [2J. Odd-dimensional 
generalizations of the formerly 4-dimensional EH metric were discovered a few years ago 
in Einstein gravity and are referred to as Eguchi-Hanson solitons. Originally obtained in 
5 dimensions upon taking a particular limit of Taub-NUT space [3], they can also be ob- 
tained from a consideration of static regular bubbles with AdS asymptotic behaviour [I] and 
generalize to any odd-dimensionality [5]. These latter cases can be derived from a set of 
inhomogeneous Einstein metrics on sphere bundles fibred over Einstein-Kahler spaces [51 E] • 

EH solitons are asymptotic to AdS^/Zp (p > 3) and have Lorentzian signature. As such 
they furnish interesting non-simply connected background manifolds for the CFT boundary 
theory. Perturbatively the EH soliton has the lowest energy in its topological class, suggesting 
that it is either a new ground state in string theory or that some as yet unknown ground 
state exists. 

Quantum gravitational corrections generally induce higher-order curvature terms in the 
low-energy effective action [8] . When Einsteinian gravity is considered in higher dimensions 
n > 4, the applicability of the theory still stands but it is not sufficient in realizing the 
free-field dynamics involved [Dl. Consequently much attention to this end has concentrated 
on Lovelock gravity [TQ1 [III E2], since its field equations are always of 2nd order, even 
though the associated action is non-linear in the curvature tensor. The most commonly 
studied special case of Lovelock gravity is Einstein-Gauss-Bonnet (EGB) gravity, which only 
contains curvature-squared terms, and implications for AdS / CFT of incorporating this term 
have recently been carried out [H]. 

In this paper we investigate the existence of EH soliton solutions in EGB gravity in 
five dimensions. We find that generalizations of the 5-dimensional EH soliton [3] do indeed 
exist. We construct such solutions via semi-analytic and numerical methods. These solutions 
smoothly approach the Einsteinian EH soliton solutions in the limit that the Gauss-Bonnet 
parameter a — > 0. 
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The outline of our paper is as follows. In Section |2j we will introduce the reader of the 
EH soliton. In Section [3j the action formalism of EGB gravity will be introduced and EGB 
field equations will be shown in Einstein-notation. Using a similar notation to [21 [2] for the 
EH soliton, we transform the metric to a more convenient form in Section |3j and also present 
the explicit EGB field equations scaled in dimensionless units. We then compute large— r 
and small-r power-series solutions to the equations in Section [4], and discuss the required 
metric regularity conditions for a soliton solution. The full numerical results of solution are 
given in Section |5j Finally, Section [6] contains a discussion of our results. 

Throughout this paper, we will use both greek and latin characters for upper and lower 
indices of the tensors but reserve latin characters for 'dummy' variables. These indices will 
run through n variables where n = 5 is the dimensionality of the spacetime. 



2 The Eguchi-Hanson Soliton 

In this section we briefly review the structure and properties of the EH soliton. Its metric 
has the form [3] 

ds 2 = -g(r)dt 2 + ^-^W + cosOdcj)] 2 + . + -dtl 2 2 (1) 

4 1 f{r)g(r) 4 2 

where d£l 2 = [dO 2 + sin 2 Od(f> 2 \ is the unit 2-sphere and 

9<r) = l + ^, /(r) = l-£, A=-| (2) 

The EH soliton is an exact solution to the 5-dimensional Einstein equations. In the I — > oo 
limit, the 5D EH soliton Q recovers the EH metric 

ds 2 = — '4^- m + cos Odcf)} 2 + — + -dQ 2 2 (3) 
4 J f(r) 4 2 v ' 

on a constant t hypersurface. 

In general the metric (II]) will not be regular unless certain conditions are imposed on its 
parameters. As r — > a, regularity in the (r, ip) section implies that ip has period In j ^Jg~{a). 
Removing string singularities at (8 = 0, it) implies that an integer multiple of this quantity 
must equal 4n. Consequently 

a 2 =( 2 {^~ l) (4) 
where p > 3 is an integer, implying that a > £. The metric ([!]) can be rewritten as 



r \ . o r 



2 



V 



tU ~ = - ( 1 + p ) dt + ~2f( r ) [ d 4> + 2 cos ( + (1 - ^ /(r) + J d ^2 (5) 



2 dr 2 r 2 



2 



where now 



f(r) = 1 - 




and we have rescaled ip so as to have period 2tt. 

The regularity condition Q implies that the metric ^ is asymptotic to AdSs/Z p where 
p > 3. The AdS/CFT correspondence conjecture states that string theory on spacetimes 
that asymptotically approach AdSs x S 5 is equivalent to a conformal field theory (CFT) 
(N = 4 super Yang-Mills U(N) gauge theory) on its boundary (S* 3 x R) x S 5 . Here the 
boundary of the EH soliton is a quotient of AdS, implying the existence of extra light states 
in the gauge theory. The density of low energy states is not affected even though the volume 
of S 3 has been reduced to S 3 /Z p . 

As noted in the introduction, the EH soliton perturbatively has the lowest energy in its 
topological class. In the context of string theory, the implications are either [5]: 

1. If the soliton is non-perturbatively stable, then it is a new ground state to string theory 

2. Otherwise, if the soliton is non-perturbatively unstable, then there exists another un- 
known ground state to string theory 

Moreover, the perturbed solutions from the EH soliton that are still asymptotic to AdSs/Z p 
(p > 3) have higher energies than the EH soliton which suggests its position in a energetically 
local minimum. 



3 Einstein-Gauss-Bonnet Gravity 



The unique combination of curvature-squared terms that is contributed by Gauss-Bonnet 
gravity towards the Einstein-Hilbert action is 



GB 



p jjabcd a p pab , p2 
tiabcdti — QKabK + K 



So that the total action Xegb becomes 

Zegb 



16ttG 



d n x^g [-2A + R + aC GB ] 



(6) 



(7) 



where n is the dimensions of the metric and a is the Gauss-Bonnet coefficient for controlling 
the relative magnitude of the contributing curvature-squared terms. If the action from ^ is 
varied with respect to the metric (and neglecting boundary terms) we obtain the EGB field 
equations 



R 



RQuv + Ao 



a 



abed 



j^abed 



AR ab R ab + R z ) - 2RR 



Liu 



+^RfiaRt + ^R^R^avb 



r\ pabc p 
Z-Tt^ J^uabc 
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Eqn ^ shows explicitly the contribution of the Gauss-Bonnet terms towards the pure Ein- 
stenian field equations. 

In order to search for EH-type soliton solutions to the equations ([8|, we employ the ansatz 



ds 2 = - r —g(r)dt 2 + - '//r • <•<>*"• 



+ 



dr 2 



+ — dfl 2 



i- 4 L " T ' TJ ' r 2 f{r)h(r) ' 4 2 

and require that {f(r),g(r), h(r)} each have an asymptotic limit of 1. In the limit a — > 0, 
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(9) 



9(r) = l + ~ 2 , h(r) = l + -, f(r) = 1 - A = - p 



(10) 



recovering the metric Q. 

The EGB field equations ^ can be written explicitly into a set of ordinary differential 
equations using We obtain 7 non-trivial field equations {£ tt , 8 rr , 8^, ^e, 8^, 8.^, 8^}. 
that in turn can be reduced to 3 independent field equations {S rr , 8^, See}- Rescaling the 
parameters and variables in the field equations via 



x 



fxx(x) = £ 2 f rr {r), f x {x) = £f r (r), A, = £ 2 A, a, 



'ID 



where the subscript denotes the derivative with respect to the relevant variable, we find 



12 fh + 12a; 4 f x h x fh - Q4fh - 32h x f + 8xhf - 128 f x h + 24x 2 f 2 h - lQx 2 f x h 



-32x 2 f xx h + 8x*f l x h l + A8x 2 f 2 h 2 + 2Ax z f xx fh + 12x z f x h x f + 80x f x fh 
+8x A f xx fh 2 + Ux 3 f x fh 2 + 24x 3 h x f 2 h\ a, 

+ 

16 - 2Ax 2 fh - 16x 3 f x h - x 4 f x h x - 6x 3 h x f - 2x 4 f xx h - 4/ - 4A*x 2 



(12) 







- A8f 2 g 2 h + 32x fh x g 2 - 2Axf 2 h x g 2 - 16x 2 fg 2 x h + 32x f x g 2 h - A8x 2 f 2 g 2 h 2 

+Ufg 2 h + Ax 4 f 2 g 2 x h 2 - 2Axff x g 2 h - 2Ax 3 ff x g 2 h 2 + \2x 2 fg 2 x h + 32x 2 jg xx gh 
-72x f 2 g x gh + 96x fg x gh - 8x A fg xx gh 2 + I6x 2 f x g x gh - A0x 3 f 2 g x gh 2 
-24x 3 f 2 h x g 2 h - \2x 2 fg x h x g + lQx 2 fg x h x g - 12x 2 ff x g x gh - 12x 4 ff x g x gh 2 



-12x 4 f 2 g x h x gh - 2Ax A f z g xx gh 



2 p2 



a* 



(13) 



4x 2 A*g 2 + 12/<f - x^fglh + 2x"fg xx gh + 2^x A fg A h + 6x 6 f x g z h 



+ 

,2 JfJ 



2f Ji 



.3/ „2i 



+6x fh x g + x f x g x gh + x fg x xh x g + lOx fg x gh - 16g 







4 



16/ 2 g 2 h - 4x 2 f 2 g 2 x h - 48x 2 f 2 g 2 h 2 + 8xf 2 h x g 2 + Ax A f 2 g 2 x h 2 + 8xff x g 2 h 

+8x 2 f 2 g xx gh + 2Axf 2 g x gh - 8x^fg xx gh 2 - 6Ax 3 ff x g 2 h 2 - 2Ax 3 f 2 h x g 2 h 
+Ax 2 f 2 g x h x g + 4:X 2 ff x g x gh - AOx 4 ff x g x gh 2 - 12x 4 f 2 g x h x gh - 8x 4 f xx fg 2 h 2 
-4x 5 f 2 g x gh 2 + 2x 5 f x g 2 x fh 2 - 4x 5 g xx f x fgh 2 - Ax b f xx g x fgh 2 

-12x 4 f x h x fg 2 h - 8x 4 f 2 g 2 h 2 - 40x 3 fg x gh 2 - 6x 5 f x g x h x fghj a* ( 14 ) 

+ 

4x 2 A x g 2 - 4fg 2 + 6x 3 fh x g 2 + 24x 2 g 2 h + 16x 3 f x h 2 h - x A fg 2 x h + 10x 3 fg x gh 
+x 4 f x h x g 2 + 2x 4 f x h x gh + x*fg x h x g + 2x 4 fg xx gh + 2x 4 f xx g 2 h = 

4 Series Solutions 

We were unable to obtain an analytic solution to the EGB field equations for the ansatz 
(J9|. We therefore turn to power-series expansions and numerical techniques to find soliton 
solutions to EGB gravity. 



4.1 Power-series solutions 

There are two particularly interesting power-series for the soliton, near infinity (the large-r 
expansion) and near r = r (the near-r expansion). While not sufficient for demonstrating 
a full soliton solution these series are useful guides as to the behaviour of the solution in 
interesting physical limits. 

For large-r we find that 

en a 4 (4cra 4 — b 4 £ 2 ) a 4 £ 2 (12aa 4 — 48a5 4 + 5b 4 £ 2 ) ( 1 

}{r > 71 + r 1 2(4a - £ 2 )r 8 15(4a - £ 2 )r w + 

£ 2 b 4 6 4 (4a6 4 - £ 2 a 4 ) a 4 £ 2 {252aa 4 - 5b 4 £ 2 + 48a6 4 ) 
9\ r ) ~ + ^2 + ~a' 2(4a-£ 2 )r 8 + 45(4a - £ 2 )r 10 + V rl ° 

1,2 b 4 b 4 (-7£ 2 a 4 + 12ab 4 ) a 4 £ 2 (756aa 4 - 5b 4 £ 2 + 336afe 4 ) / 1 



lrj + r 2 + r 4 6(4a - £ 2 )r 8 45(4a - £ 2 )r w V 10 

(15) 

where 4a - £ V (4a, -1^0). 

There are two free parameters (a 4 , 6 4 ) in the solution, governing the falloff rate of the 
metric functions. Their falloff rates suggest that the mass of the soliton, M, is governed by 
these two quantities. Using the conformal formalism based on the electric Weyl tensor of 
[T3| [14] to compute the conserved quantities, the conformal mass in Einstein gravity is given 
by 

7r(36 4 — a 4 ) 



8£ 2 Gp 



(16) 
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Setting a 4 = — r 4 and 64 = in (16), we recover the soliton mass from the counterterm 
subtraction method [3 J and from the Hamiltonian formalism [I] . 

To obtain the near-r power-series solution, we impose the conditions /(r ) = 0, g(r ) 7^ 0, 
and h(ro) 7^ at the bubble-edge tq. We then find 

f(r) = A x {r - r ) + A 2 (r - r ) 2 + A 3 (r - r ) 3 + O ((r - r ) 4 ) 

g{r) = B + B x {r - r Q ) + B 2 (r - r ) 2 + B 3 (r - r ) 3 + O ((r - r ) 4 ) (17) 

/i(r) = Co + C 1 (r-ro) + C 2 (r-r ) 2 + C 3 (r-ro) 3 + 0((r-ro) 4 ) 
6 12a 

where A = — — H — — - and B 7^ and Co 7^ 0. The coefficients in the preceding expression 
are extremely cumbersome and we relegate them to Appendix [B} 

Before moving to a numerical solution of the equations, we revisit the issue of metric 
regularity for the EGB case. The Kretschmann scalar /C = R abca ''R abcdl 

/C = — ^ (Vr 8 /^ 2 - 4r 8 gh 2 fg rr g 2 r + 32 g 4 r 7 f r h 2 f rr + Ar 8 g 2 h 2 fg rr g r f r 

+4r 8 g 2 hf 2 g rr g r h r + lQgWff rr + 16rV^/ 2 + r 8 /* 2 A 4 + 88rV#i a 
+ 160rV Z^ 2 / 2 + If* W/V + 8r 7 g 3 h 2 f g„ f r + Ar 8 g 2 h 2 f 2 g 2 „ 
+2Ar 7 g 2 h 2 f 2 g r g rr + 8r 7 g 3 f 2 h 2 g r + Ar 7 g 3 h 2 r f 2 g r + 8r 7 g 3 hf 2 g rr h r 
+8g 4 r 7 f rr hh r f + r s g 2 g 2 h 2 f 2 + Wg 4 r 7 f 2 hh r + A8f 2 h£ 2 g 4 r 4 - 12r 7 gh 2 f 2 g 3 
+8r 7 g 3 h r fg r f r h + 2r 8 g 2 g 2 h r f f r h + Ag 4 r 8 f rr hh r f r + Ag 4 r 7 h 2 ff r (18) 
-I28g 4 £ 2 r 2 fh + 32g 4 r 2 f 2 h£ 2 + 2r 8 g 2 g 2 f 2 h 2 + A0r 6 g 2 h 2 f 2 g 2 + 64r 5 g 4 hf 2 h T 
+lQ0r 5 g 4 h 2 ff r + 96r 5 g 3 h 2 f 2 g r + 17Qg 4 £ 4 f 2 - 384g 4 f£ 4 + g 4 r 8 h 2 J 2 
+256£ 4 g 4 - 2r 8 hf 2 g 3 .gh r + 8r 7 g 2 hf 2 g 2 h r + 12r 7 g 2 h 2 f g 2 f r 
+6Ar 6 g 4 h r ff r h + A8r 6 g 3 h 2 fg r f r - 32g 4 r 3 f r hf£ 2 - 2r 8 h 2 fg 3 gf r 
+32r 6 g 3 hf 2 g r h r ) 

and will be finite for all r > r provided g(r) 7^ for all r > r . The only remaining possible 
singularities will be conical singularities at r = ro and string singularities. The former will 
not be present provided ip has a period of 27r/V where 

= r 4 A 2 R 2 C Q 

- we 1 ' 

String-singularities will be eliminated at the north and south poles if ip has a period of 4ir/p 
where p 6 Z \ {0}. Together these conditions imply 

- ■ M 
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This is a generalization of the regularity condition for the Einsteinian case [5]. For a = it 

4 f 
is straightforward to show that A\ = — and Co = 1 H — ~ giving 

r r^; 

^ = Z 2 ( J " l) (21) 
which is the regularity condition (|4]), with a = ro. 

5 Numerical Results 

In this section we present the numerical solutions of the EGB field equations. The power- 
series solutions from Section [4] provide a useful guide for approximating the initial conditions 
(hereafter IC) of the large-r and near-ro field equations. Both positive and negative values of 
the cosmological constant are possible, depending on the choice of Gauss-Bonnet coefficient 
a*. 

We have compared the full numerical solutions with their respective power-series expan- 
sion as a cross-check on our numerical work. As shown more fully in Appendix [AJ the large-r 
numerical solution agrees very well with the corresponding power series solution when ex- 
panded to order 1/r 10 , even for reasonably small r. The near-ro power-series agrees well 
with the numerical solution when the corresponding power series solution is expanded up 
to (r — ro) 3 order, but can quickly deviate from the numerical solution when r is appre- 
ciably larger than r . We use the power series expansions given in Section [4] for the initial 
conditions. 

With consistent results from the power-series approximations, we now present the full 
numerical solutions below. 



5.1 Initial Conditions and Numerical Procedure 



We find that there is numerical instability if the values of |a*| are too large. From the series 



solutions (15) we see that a* ^ 1/4, so we shall consider only small a* < 1/4 henceforth. We 



can choose freely the parameter p, which is a positive integer from condition (20) required 



to satisfy metric regularity conditions. In Einstein gravity, we have the bound p > 3 from 



the regularity condition (21). The analogous lower bound on p in EGB gravity follows from 
condition (20) since the coefficients A\ and Cq will depend on ro and i, which we do not 
know a priori. For sufficiently small a* we expect that p = 5 will be sufficient, and we shall 
make this choice throughout the paper. It is a straightforward (albeit tedious) exercise to 
numerically obtain solutions for larger p. 

All other parameters are constrained by the shooting method, a trial-and-error method 
of tuning the values of chosen IC parameters that are not known a priori, so that when the 
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equations are numerically integrated away from the starting point, the solution will satisfy 
desired asymptotic properties. The IC parameters {a* A , b\} and {B , C } will be used to 
tune the required desired properties of the solution for the large-r and near-ro solutions, 
respectively. For numerical integration beginning near-ro, we use the shooting method to 
require {f,g, h} to follow the required large-r asymptotic conditions. Conversely, for large-r 
numerical integration, the starting point is at r = oo, and the shooting method is used to 
set the value of the bubble edge at r = r . 

For convenience, we list here the various parameters we employ for both large-r and 
near-r numerical integration of the field equations. 



• a* 

This is the dimensionless Gauss-Bonnet coefficient, chosen to be small to avoid issues 
of numerical instability as noted above. 



P 

This is a non-zero integer that emerges from the the metric regularity condition (20), 
whose choice selects one of a countably infinite set of EH soliton solutions in EGB 
gravity. The 5D Einsteinian EH soliton required that p > 3 [3] . For practical purposes, 
we will use p = 5 for our numerical results. 

Bo 

This the value of g{r = r ); it must be positive and converge to 1 + ^ 2 / r o in the a — > 
limit. The near-r shooting method will constrain this parameter. 

C 

Similar to B , this is exactly the value of h(r = r ) and should be positive and con- 
verging to 1 + £ 2 / r\ in the a — >■ limit. The near-ro shooting method will constrain 
this parameter. 



• Xo 

This is the dimensionless radius of the EH soliton bubble edge ro in units of I defined 
by xo = tq/£. The spacetime exists only for x > xq (equivalently r > ro). The large-r 
shooting method will constrain this parameter. 



This parameter contributes to the mass of the soliton. In the Einsteinian limit a — > 0, 
(X4 = — r , where r is the soliton radius. This interpretation will not necessarily hold 
for nonzero a and 64. The large-r shooting method will constrain this parameter, whose 
rescaled version is a* 4 = a^ji . 

h 

The 64 parameter also contributes to the mass of the soliton; in Einstein gravity the 



relevant expression is that of Eqn (16). We will consider both 64 = and 64 7^ cases. 
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The large-r shooting method will constrain this parameter, whose rescaled version is 
K = h/£ A . 

• so 

This is the reciprocal dimensionless radius of the EH soliton bubble edge, defined by 
s = £/r . To numerically analyze the r — > oo limit, we make a variable change 
r (->■ s — £/ro and analyze its s — > behavior. Note that the spacetime is valid for 
s < so only. 

The outline of the numerical procedure is as follows. For some a* small, we choose p — 5. 
We then make an initial choice of the IC parameters {a* 4 , b* 4 } (large-r) or {B , Co} (near- 
r ). The other parameters in the respective power-series solutions are then fully determined 
from this guess. We then apply the Runge-Kutta-Fehlberg method to integrate the field 
equations £ rr ,£^,£gg, with the starting point at r = r (near-r ) or at r = oo (large-r). 
If the required asymptotic behaviour of the soliton (when starting at r = r ) or bubble 
edge behaviour (when starting at r = oo) is satisfied, we are done and have obtained the 
numerical solution. If they are not satisfied within a given tolerance (we use a tolerance of 
e = 10~ 4 ), then we repeat the process with improved guesses of the ICs and calculate the 
numerical solution. 

5.2 Large-r behavior 

We employ the parameter s = ijr < 1 to investigate the large-r behavior of our solutions. 
Our choice of p = 5 implies that Sq =0 = l/x9. = (p 2 /4: — l) 0.4364 in the Einsteinian 

p=5 

case. In Einstein gravity, we also have a* 4 ' a ~° = —x^ = —27.5625 and b* A ,a ~° = 0. We find 

p=5 

that the field equations for large-r are numerically more difficult than near-r , so we only 
display results that we numerically trust. 

The free parameters chosen for our large-r numerical results are given in Table 1. We 
present our results for the parameter set fi G {a* 4 , b* 4 , Sq} in terms of their deviations A/i = 
/i — fi a=0 relative their values in Einstein gravity. 
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Numerical Parameters — large-r 


a* = 


a* = -0.04 


a* = 0.002 


Figure Number 




#2 


A3 






P 


Free 


5 


5 


5 


5 


5 


Aa* 


Shooting 
method 





7.80 


7.80 


-0.20 


-6.30 


A6: 


Shooting 
method 








-2.96 





0.200 


As 


Calculated 





5.17 x 10~ 2 


2.00 x 10~ 2 


2.20 x 10~ 3 


-2.11 x 10~ 2 


AM* 


Calculated 





1.56 


3.34 


-0.04 


-1.38 



Table 1: Numerical values for the large-r numerical solution. The symbols used are sq 



1/xq, s% =0 « 0.436436, a 4 



,o=0 



-27.5625, 



*,a=0 



the a = values are defined by Aa^ 
The mass deviation parameter AM* 



0, and 
Abl 



*' a = —5.51. Deviations from 

*,a=0 



" !l *bl = b\ - 6r=°, and As = s - s%=°. 
p,a=o j g j n un j^ s Q f Til 2 /%Q^ where the latter 



term is given in Eqn (16). 



We find that the numerical integration is much more sensitive to the choice of h\ than a\, 
and so require more fine-tuning of b\. 

Figure [I] - Figure [5] depict the solutions {f(s),g(s),h(s)} where s = ijr for a* = 
0, —0.04, 0.002 from Table [TJ respectively. We see that all metric functions approach 1 
in the s — > limit (at r — > oo) as required. Our solutions are valid all the way from infinity 
to the edge of the soliton, whose values are given in Table [T] for the various cases. 

We also see from Fig. [3] that a nonzero value for 64 can have a dramatic effect relative 
to the 64 = 0, a = —0.04 solutions shown in Fig. [2] This effect is less pronounced for the 
a = 0.002 solutions, as shown in Fig. |4]and Fig. [5]possibly due to the decreasing magnitude 
of a. 
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2 -i 



1.5 - 



0.5 - 








£ 

Figure 1: Large-r numerical solutions {/, g, h} for a* = where s = -. We have a\ = a* 4 ' a 

_ r 
64 = and so = Sq =0 here. The color coding is f(s) = blue, g(s) = red, and h(s) = green. 




s 



Figure 2: Large-r numerical solutions {/, g, h} for a* = —0.04 where s = -. We have 

r 

a\ = a* 4 ,a ~° + 7.80, b\ = 0, s = Sq =0 + 0.0517 here. The color coding is f(s) = blue, 
g(s) = red, and h(s) = green. 
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2 -i 



1.5 - 



1 - 



0.5 - 








s 

£ 

Figure 3: Large-r numerical solutions {f,g,h} for a* = —0.04 where s = -. We have 

_ r 
a \ = al' a=0 + 7.80, b\ = -2.96, s = s% =0 + 0.02 here. The color coding is f(s) = blue, 

g(s) = red, and h(s) = green. 

2 -i 



1.5 - 



1 - 



0.5 - 




o H 1 1 1 1 1 1 ■ r- 

0.1 0.2 0.3 0.4 



s 

Figure 4: Large-r numerical solutions {f,g,h} for a* = 0.002 where s — -. We have 

_ r 

a* = al' a=0 - 0.20, b\ = 0, s = s% =0 - 0.0022 here. The color coding is f(s) = blue, 

g(s) = red, and h(s) = green. 
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Figure 5: Large-r numerical solutions {/, g, h} for a* = 0.002 where s 



We have 



6.30, b\ 



9[s 



= 0.20, s 
green. 



0.0211 here. The color coding is f(s) = blue, 



red, and h(s 
5.3 Near-ro behavior 

The free and calculated parameters with their numerical values for the near-ro solution 
are shown in Table § In the Einsteinian case, B* =0 = C Q=0 = 1 + (l/< =0 ) 2 = 1 + 

(p 2 /4-l)~ 1 « 1.1905 where x% =0 = (p 2 /4-l) 1/2 « 2.2913. As in the large-r case, we 

p=5 p=5 

express our results in terms of deviations from Einstein gravity. 





Numerical Parameters — near-r 


a* = 


a, = -0.04 


a, = +0.002 


Figure Number 


fl6 


fl7 




V 


Free Parameter 


5 


5 


5 


AB 


Shooting method 





-0.1423 


0.0087 


AC 


Shooting method 





-0.1782 


0.0065 


Ax 


Calculated 





-0.2442 


0.0116 



Table 2: Numerical values for the near-ro solution, where B t 



a=0 



o=0 




Ga=0 




1.1905 and 



= 2.2913. 

C 



Deviations from the a = values are defined by ABq = Bq — B s 



a=0 



C^ , and Ax c 



x 



x, 



a=0 
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The only free parameter in the near-r solution is p; all other parameters can be obtained 
numerically via the shooting method and from solving the constraint equations. 
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Figure 6: Near-r numerical solutions of {/, g, h} for a* = 0, p = 5, B 



and x 



x, 



a=0 



r>a=0 
B ) 



Co 



Ca=0 
) 



The color coding is f(x) = blue, g(x) = red, and h(x) = green. 
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Figure 7: Near-r numerical solutions {f,g, h} for a* = —0.04, p = 5, B 



C 
and h(x 



°0 



- 0.1782, and x 
green. 



x 



a=0 



0.1423, 



0.2442. The color coding is fix) = blue, g(x) = red, 
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0. 




"I— 1 ■ 1 ■ 1 1 1 ■ 1 r- 

10 20 30 40 

x 

Figure 8: Near-r numerical solutions {/, g, h} for a* = +0.002, p = 5, B = Bq =0 + 0.0087, 
C = C£ =0 + 0.0065, and x = x% =0 + 0.0116. The color coding is f(x) = blue, g(x) = red, 
and h(x) = green. 

6 Conclusions 

We have semi-analytically and numerically found five- dimensional Eguchi-Hanson soliton 
solutions in Einstein-Gauss-Bonnet (EGB) gravity. We have illustrated the typical form of 
the metric functions for small positive and negative values of a, integrating from both large- 
r to the edge of the soliton, and from the edge of the soliton to infinity. These numerical 
solutions are fully consistent with large-r and the near-r power-series solutions with 2 free 
parameters a and p. 

We have found numerical evidence (in the context of the large-r problem) that 64 can be 
non-vanishing in EGB gravity. This indicates a broader class of soliton solutions, in which 
both 04 and 64 can be nonzero. A full exploration of these solutions, along with a study of 
a broad range of a* (going beyond the choice we believe falls in a safe range of stability), 
remain interesting subjects for future investigation. 
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7 Appendix 



A Numerical Solutions Compared to Power Series Ex- 
pansions 

We present here a comparison of our numerical solutions to the power-series expansion of 
our results, up to 1/r 10 for the large-r solution and up to (r — r ) 3 in the near-r solution 
for a = 0. 

The large-r power-series is very effective when expanded up to 1/r 10 , but the near-r 
solutions fail to demonstrate proper asymptotic behavior when expanded only to order (r — 
ro) 3 . 
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(a) The large-r power-series for a* = 0. 



(b) The large-r numerical solution for a* = 0. 




(c) The large-r power-series for a* = —0.04. (d) The large-r numerical solution for a* = —0.04. 



Figure 9: The large-r solution comparison of the power-series in expanded up to 1/r 10 
(dashed lines) and full numerical solution (solid lines) for a* = and a* = —0.04 with a 
colour coding of / = blue, g = red, and h = green. 
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Figure 10: The near-r solution comparison of the power-series in expanded up to (r — r ) 3 
(dashed lines) and full numerical solution (solid lines) for a* = and a* = —0.04 with a 
colour coding of / = blue, g = red, and h = green. 
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0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 

x x 

(a) The deviation for large-r solutions for a* = 0. (b) The deviation of the large-r solutions for a* = 

-0.04. 




(c) The deviation of the near-ro solutions for a* = 0. (d) The deviation of the near-ro solutions for a* = 

-0.04. 

Figure 11: The deviation, A = F ps (x) — F num (x), of the large-r and near-ro power-series 
solution with the full numerical solution for a* = and a* = —0.04 with a colour coding of 
/ = blue, g = red, and h = green. The large-r power-series is expanded up to 1/r 10 , and 
the near-r power-series is expanded up to (r — r ) 3 . 
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B Large-r power-series coefficients 



We present here the near-r power series coefficients up to (r — r ) from Eq. (17) 



Ui 2 P 2 



A, 



A* 



An 



A? 
A? 



5 r w tpC + (2£ 5 y^Co - 68 £ 3 p 2 C 3/2 a) r 9 
+ (-96 f P C a - fp + 336 £ 2 p 3 C 2 a 2 ) r 8 

+ (-704£p 4 C 5/ V - 20 £ 5 p 2 a ^/C~ + 672 £ 3 p 2 C 3/2 a 2 - 16£ 5 a ^/C^j r 7 

+ (-24 fpa + 512aVC 3 - 1280 £ 2 p 3 C 2 a 3 + 224 fp 3 C a 2 - 7G8 £ 4 pC a 2 ) r 6 

+ (-768 a 2 £ 5 y/C~ + 8448 £ 3 p 2 C 3/2 a 3 + 736 ^pV ^/C~ - 768 « 3 ^ 4 C 3/2 ) r 5 

+ (l024a 4 £VC 2 + 12288 £ 4 pC a 3 - 20480 « 4 lVC 2 - 5 760 a 3 £ A p 3 C + 16 a 2 fp 3 ) r ' 

+ (-2816£Va 3 v / C^+ 14336a 4 £VCo 3/2 - 36864 a 4 fp 2 C 3/2 - 64a 3 £V v 7 ^) r 3 

+ (14336 aH A p 3 C - 1536a 4 ^yC - 1408 a 3 £ 6 p 3 + 2048 a 3 £ 6 p) r 2 

+ (2048 a 4 £V - 20480 aH b p 2 ^/C^j r + 10240 a 4 £ 6 p 3 



r 5 C 3/2 (r 2 + 16 a) (4 r apy/&) + 8 a£ - £ r/) (64a 2 £p v / ^) r o 
+16p 2 C a 2 r 2 - 16a 2 £V - 16r 2 £ 2 a - 8£ r 3 ap^[C^ + r 4 £ 2 ) 



= -AB £ 7 ^ r 13 + 64:B a£ 6 pC ro 12 -2AB £ 5 a >/C~ (-4£ 2 + £ 2 p 2 + 8p 2 C a) r x 
-32 B a 2 £ 4 pC (-15£ 2 p 2 + 12£ 2 + 80p 2 C a) r 10 

+64 B a 2 £ 3 ^ (17 £ 4 p 2 + 400p 4 C V + 8£ 4 - 240 £ 2 p 2 C a - 60aiYC ) r 9 
-1024fi a 3 ^ 2 pCo (-15a£VC + 34 1 + 96p 4 C V - 180£ 2 p 2 C a + lQfp 2 ) r 8 
+25QB a 3 ^/C~ £ (704 a 3 p 6 C 3 + 185Qp 2 a£ 4 C -9Qf- 52 tp 2 - 120 a 2 £ 2 p 6 C 2 
-3360 a 2 £ 2 p*C 2 + 344 a £ 4 p 4 C + 3 £ 6 p 4 )r 7 
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-1024 B a 4 pC (17Q a £ 4 p 4 C - 24 a 2 £ 2 p 6 C 2 + 9 tp 4 + 128 a 3 p 6 C 3 - 480 1 
-182A a 2 £ 2 p 4 C 2 - 124 fip 2 + 2624 p 2 a £ 4 C )r 6 

-1024 B a 4 £ ^/C~ (-128£ 6 + 240a^p 4 C + 1536 «VC 3 + 3456p 2 a£ 4 C - 16£ 6 p 2 
-6784« 2 £VC 2 - 36a^p 6 C - 16« 2 £VC 2 + 19 fp 4 )r Fj 
+lQ384B a 5 £ 2 pC (12£ 4 p 2 + lQC 2 p 6 a 2 -9Q£ 4 - 41Qp 4 C 2 a 2 + 672 £ 2 p 2 C a 
-3C p e £ 2 a + 9p 4 £ 4 - Ma£ 2 p 4 C )r 4 + 2048 B a 5 £ 3 p 2 ^(-1152 £ 2 p 2 C a 
+U£ 4 p 2 - 61Up 2 C 2 a 2 - 3p 4 £ 4 + 3072 £ 2 C a + 1024p 4 C V + 128 f 
-120 a £ 2 p 4 C )r 3 

-8192 B a 6 £ 4 p 3 C (Sfp 4 + 256 £ 2 + 28£ 2 p 2 + lQap 4 C - Q40p 2 C a + 1024 a C ) r 
-16384 B a 6 £ 5 p 4 ^/C~ (-256 aC + 32p 2 C a - 3 £ 2 p 2 - 8 £ 2 ) r - 524288 B a 7 p 5 C £ 6 

r p (r 2 + 16 a) (a r a py^C~ + 8a£-£ r 2 ) 2 (64 a 2 £p^/C~ r + 16 p 2 C a 2 r 2 
-16 a 2 £ 2 p 2 - 16 r 2 £ 2 a -8£ r 3 ap^ + r 4 £ 2 ) 2 C 



+8a£ 2 >/C~ (l5£ 4 p 2 - 248 £ 2 p 2 C a - lUa£ 2 p 4 C + A£ 4 + 832p 4 C V) B r 9 
-4a£p(20f - 944a 2 £VC 2 + 3£ 6 p 2 - 3392p 2 £ 2 C V + 640 £ 4 C a 
+5Q0p 2 a£ 4 C + 48Q4p 4 a 3 C 3 )B r 8 



+32a 2 v ^C~ (2A0p 2 a£ 4 C - 576 a 2 £ 2 p 4 C 2 - 160a 2 £VC 2 + 340 a £ 4 p 4 C - 64 £ 6 
+640 a 3 p 6 C 3 + 8 fp 2 + fp 4 )B r 7 

-256a 2 £p{320p 4 a 3 C 3 + fp 2 - 592 p 2 £ 2 C 2 a 2 + 3£ 6 + 30 a 2 £ 2 p 4 C 2 



+7Ap 2 a£ 4 C + A8£ 4 C a - 2 a£ 4 p 4 C )B r 6 + 128 a 3 ^(-140^ + 43 
-192£ 6 - 14a£VC + 3456p 2 a£ 4 C + U16a£ 4 p 4 C - 224a 2 £VC 2 
+ 1536 a 3 p 6 C 3 - 7552 a 2 £ 2 p 4 C 2 )B r 5 

+Q4a 3 £p(-8Q4a£ 4 p 4 C - 37888 p 2 £ 2 C 2 a 2 - 36 tp 2 + 4352p 2 a£ 4 C 
+6144£ 4 C a - 5760 a 2 £ 2 p 4 C 2 + 64 £ 6 + 3£ 6 p 4 + 24576 p 4 a 3 C 3 )B r 4 



-25Qa 4 £ 2 p 2 ^ (-15360p 2 C 2 a 2 + 236 l 4 p 2 + 2560 /VC a + 7680£ 2 C a 
+21p 4 £ 4 + 768p 4 C V - 128 1 - 560 a £ 2 p 4 C )B r 3 
-512a 4 £ 3 p{768£ 2 p 2 C a - 128 £ 4 + 1536p 4 C V - 512a^VC + 33p 4 £ 4 
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-61Up 2 C 2 a 2 )B r 2 

-4096 a 5 £YVCo (l60£ 2 - 28 £ 2 p 2 + !92p 2 C a - 15 £ 2 p 4 ) B r 
+40960 aH 7 p 3 (3p 2 + 8) B 

= r 3 C 3/2 p 2 (Q4a 2 £p^/Q ) ro + 16p 2 C a 2 r 2 - 16 a 2 £ 2 p 2 - 16 r 2 £ 2 a 
-8£r 3 apy^C~o + r 4 £ 2 ) (4 r apy/C~ + 8 al - £ tq 2 ^ (r 2 + 16a) 2 



Cf = 2 £ 3 p (f + 16 a C ) r 8 - 16 a £ 2 y^C~ (-2 £ 2 + 24 p 2 C a-£ 2 p 2 )r 7 
+ 16a£p(-16£ 2 C a- 14 £ 2 p 2 C a + 3 £ 4 + 96 p 2 C 2 a 2 ) r 6 

-256 a 2 t/C~ (-2 1 + 8 p 4 C V - 2 a £ 2 p 4 C + 3 £ 4 p 2 - 2 £ V C a) r 5 
+32 a 2 £ 3 p (-£ 2 p 2 + 104p 2 Q,a - 128 a C ) r 4 

-512 a 3 £ 2 p 2 y^C~ (-16 a C + 8 p 2 C a + 3 £ 2 -£ 2 p 2 ) r 3 
+256 o?£ 3 p (ll£ 2 p 2 - 16£ 2 + 48p 2 C a) r 2 

-10240 a 4 £ 4 p 2 Vc~o(-2 + p) (p + 2) r - 20480 a 4 £ 5 p 3 

Cf = r 3 p(6Aa 2 £p^Cor + 16p 2 C a 2 r 2 - 16a 2 £V - 16 r 2 £ 2 a 

-8£r 3 ap^Co + r i £ 2 ) (a r ap^/Co + 8 a £ - £ r 2 ) (r 2 + 16 a) 



cl 
c» 

= -2£ 9 r 22 pC + 2£ 8 ^C~ (l44p 2 C a + 2£ 2 + 5£ 2 p 2 ) r 21 
-2£ 7 p (8£ 2 C a + £ 4 + 4032 p 2 C 2 a 2 - 54£ 2 p 2 C a) r 20 

+8a£ 6 ^/C~ (-942 a £ 2 p 4 C + 13440 p 4 C 2 a 2 - 328 £ 2 p 2 C a - 20 f + 71£ 4 p 2 ) r 
-Aa£ 5 p{-39fp 2 + 209664 p 4 a 3 C 3 + 6456p 2 a£ 4 C - 14272 p 2 £ 2 C 2 a 2 
-2A6A£ 4 C Q a - 4£ 6 - 27216 a 2 £ 2 p 4 C 2 )r 18 

+ 16a 2 £ 4 v / C^(258048aVC 3 - 96 £ 6 - 292 ^p 2 - 34496 a 2 £ 2 p 4 C 2 - 223 fp 4 
-49392 a 2 £ 2 p 6 C 2 + 27720 a £ 4 p 4 C - U112p 2 a£ 4 C )r 17 
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-64a 2 £ 3 p(-716p 4 £ 6 aC + 61704 p 4 1 a 2 ] C 2 - 2280p 2 £ 6 aC - 52080 p 6 £ 2 a 3 C 3 

-llp 2 £ 8 - 288£ 6 aC + 204288 p 6 a 4 C 4 - 4QU4p 4 £ 2 a 3 C 3 - 32£ 8 

-41504 p 2 £ 4 a 2 C 2 )r 16 + 128 a 3 £ 2 ^(-65184 a 3 £ 2 p 8 C 3 + 704 £ 8 

+ 160176 a 2 £ 4 p 6 C 2 - 12528p¥ 5 aC - 183Qp 2 £ 8 - 4352p 2 faC + 439 £ 8 p 4 

-Q8992p 6 £ 2 a 3 C 3 - 32Q4a£ 6 p e C - 140992 p 4 £ 4 a 2 C 2 + 20 2 752 a 4 p 8 C 4 )r 15 

-Q4a 3 £p(-lU4320p 4 £ 4 a 3 C 3 - 211968 p^VQ, 2 - 185 3 44 p 6 £ 2 a 4 C 4 

-18201Q p 8 £ 2 a 4 C 4 + 10087Q8p 6 £ 4 a 3 C 3 + 51456 £ 8 aC - 123200 p 4 £ 6 a 2 C 2 

+460800 p 8 a 5 C 5 - 39920p e fa 2 C 2 + 660 p 2 £ 10 - 21p 4 £ 10 + 29320 p 4 £ 8 aC 

-9385Qp 2 £ 8 a C )r 14 + 256 a 4 y / C~ (-275904p 4 £ 8 a C + 479616 a 3 £ 4 p 8 C 3 + 512 £ 10 

-37744 a 2 fp 8 C 2 + 20 7360 p 2 £ 8 a C - 7772 p 4 £ 10 + 57344 a 5 p 10 C 5 

-591360 p 4 £ 6 a 2 C 2 - 28416 a 4 £ 2 p w C 4 + 56 00 p 2 £ 10 - 70464p 6 £ 6 a 2 C 2 

+85328 afp 6 C - 669184p 6 £VQ, 3 - 59£ 10 p 6 + 22528 p 8 £ 2 a 4 C 4 )r 13 

+ 1024a 4 £p(-3Q8Q4p 8 a 5 C 5 + 20 4 48 a 3 £ 4 p 8 C 3 + 49 8 04 p¥a(7 

- 134400 p 8 £ 2 a 4 C 4 - 134 p 4 £ 10 + 33792 p 6 £ 4 a 3 C 3 + 866816 p 4 £ 4 a 3 C 3 
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-389120p 6 £ 2 a 3 C 3 + 1894912 p 4 £ 4 a 2 C 2 + 24 5 76 p 2 £ 4 a 2 C 2 + 30780 £ 8 p 4 

+ 1004800 a 2 £ 4 p 6 C 2 + 44 2 3 68 £ 6 aC + 48640 £ 8 + 851968 p 6 a 4 C 4 - 10656 p 2 £ 8 

- 1641984 p 2 fa C ) r 9 - 16384a 6 £ 3 p(-320p w £ 2 a 3 C 3 + 4096p 10 a 4 C 4 - 

1766400 p 6 £ 2 a 3 C 3 - 290816 p 2 £ 6 a C + 2435072 p 4 £ 4 a 2 C 2 - 20 48 £ 8 
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-32768 a 7 £ 4 p 2 v / C~ (32768p 4 a 3 C 3 + 393216 p 2 £ 2 C 2 a 2 - 988 a £ 4 p 8 C 
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References 

[1] J. Maldacena. The Large-N Limit of Superconformal Field Theories and Supergravity. 
International Journal of Theoretical Physics, 38:1113-1133, 1999. 

[2] A. M. Ghezelbash and Robert B. Mann. Atiyah-Hitchin M-branes. JEEP, 10:012, 2004. 

[3] R. Clarkson and Robert B. Mann. Soliton solutions to the Einstein equations in five 
dimensions. Phys. Rev. Lett, 96:051104, 2006. 

[4] Keith Copsey. Bubbles Unbound II: AdS and the Single Bubble. JEEP, 10:095, 2007. 

[5] R. Clarkson and R. B. Mann. Eguchi Hanson solitons in odd dimensions. Classical and 
Quantum Gravity, 23:1507-1523, March 2006. 



24 



D. N. Page and C. N. Pope. Inhomogeneous Einstein metrics on complex line bundles. 

Classical and Quantum Gravity, 4:213-225, 1987. 

H. Lii, D. N. Page, and C. N. Pope. New inhomogeneous Einstein metrics on sphere 
bundles over Einstein-Kahler manifolds. Physics Letters B, 593:218-226, July 2004. 

N. Ogawa and T. Takayanagi. Higher Derivative Corrections to Holographic Entangle- 
ment Entropy for AdS Solitons. ArXiv e-prints, July 2011. 

N. Dadhich. On the Gauss-Bonnet Gravity. ArXiv High Energy Physics - Theory 
e-prints, September 2005. 

D. Lovelock. The Einstein tensor and its generalizations. J. Math. Phys., 12(498), 1971. 

D. Lovelock. Divergence-free tensorial concomitants. Aequantiones Math., 4(127), 1970. 

D. Lovelock. Tensors differential forms and variational principles. Wiley-Interscience, 
New York, 1975. 

A. Ashtekar and S. Das. LETTER TO THE EDITOR: Asymptotically anti-de Sitter 
spacetimes: conserved quantities. Classical and Quantum Gravity, 17:L17-L30, January 
2000. 

S. Das and R. B. Mann. Conserved quantities in Kerr-anti-de Sitter spacetimes in 
various dimensions. Journal of High Energy Physics, 8:33, August 2000. 



25 



